###########################################################################
# Plot
# Figure 3

par(mar=c(4,3,2,1), mfrow=c(1,1), mgp=c(2, .5, 0), cex.lab=.8, cex.axis=.8)

plot(av.ACME.mi, 4, cex=1.5, axes=F,  ylab="", xlim=c(-3,3), ylim=c(.5, 4.5), main="", pch=17, xlab="ACME of Permitting a Demonstration on Feeling Thermometer Scores")
grid(lty=1, lwd=.5, col="lightgrey")
abline(v=0, col="grey")
axis(1, col="white")
axis(3, col="white")

points(av.ACME.ci, 3, pch=17, cex=1.5)
points(av.ACME.mn, 2, pch=17, cex=1.5)
points(av.ACME.cn, 1, pch=17, cex=1.5)
segments(av.ACME.mi-1.96*av.ACME.se.mi, 4, av.ACME.mi+1.96*av.ACME.se.mi, 4)
segments(av.ACME.ci-1.96*av.ACME.se.ci, 3, av.ACME.ci+1.96*av.ACME.se.ci, 3)
segments(av.ACME.mn-1.96*av.ACME.se.mn, 2, av.ACME.mn +1.96*av.ACME.se.mn, 2)
segments(av.ACME.cn-1.96*av.ACME.se.cn, 1, av.ACME.cn+1.96*av.ACME.se.cn, 1)

text(-3, 4.2, "Muslim Immigrants", pos=4, cex=.8)
text(-3, 3.2, "Christian Immigrants", pos=4, cex=.8)
text(-3, 2.2, "Muslim Natives", pos=4, cex=.8)
text(-3, 1.2, "Christian Natives", pos=4, cex=.8)